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ABSTRACT 

Aims. Recent theoretical predictions of the lowest very high energy (VHE) luminosity of SN 1006 are only a factor 5 below the previously pub- 
lished H.E.S.S. upper limit, thus motivating further in-depth observations of this source. 

Methods. Deep observations at VHE energies (above 100 GeV) were carried out with the High Energy Stereoscopic System (H.E.S.S.) of 
Cherenkov Telescopes from 2003 to 2008. More than 100 hours of data have been collected and subjected to an improved analysis procedure. 
Results. Observations resulted in the detection of VHE y-rays from SN 1006. The measured y-ray spectrum is compatible with a power-law, the 
flux is of the order of 1% of that detected from the Crab Nebula, and is thus consistent with the previously established H.E.S.S. upper limit. The 
source exhibits a bipolar morphology, which is strongly correlated with non-thermal X-rays. 

Conclusions. Because the thickness of the VHE-shell is compatible with emission from a thin rim, particle acceleration in shock waves is likely 
to be the origin of the y-ray signal. The measured flux level can be accounted for by inverse Compton emission, but a mixed scenario that in- 
cludes leptonic and hadronic components and takes into account the ambient matter density inferred from observations also leads to a satisfactory 
description of the multi-wavelength spectrum. 

Key words, y-rays: observations - SNR: individual (SN 1006, G327.6+14.6) - supernova remnants 



1. Introduction 

The source SN 1006 is the remnant of one of the few historical super- 
novae. It appeared in the southern sky on 1006 May 1 and was rec orded 
by Chinese and Arab astronomers IStephenson & Green 2003l . The 
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remnant of this explosion was first identified at radio wavelengths on 
the basis of historical evidence IGardner & Milne 1965H . The evolution 
of its luminosit y indicates that it is the result of a Type la supernova 
ISchaefer 19 961. probably the brightest supernova in recorded history. 
A distance of 2.2 kpc was derived by Winkler et al. (2003) based on 
comparing the optical proper motion with an estimate of the shock ve- 
locity derived from optical thermal line broadening assuming a high 
Mach number single- fluid shock. 
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Contemporary interest in the very high energy (VHE) emission 
from supernova remnants (SNRs) has arisen due to their association as 
prime candidates for Galactic cosmic-ray acceleration. Firstly, Galactic 
SNRs have sufficient kinetic energy to explain the estimated Galactic 
luminosity in cosmic rays of 10 40 erg/s. Secondly, and more impor- 
tantly, it has been shown that diffusive shock acceleration provides a 
viable mechanism which can efficiently accelerate charged particles in 
the blast waves of SNRs (e.g. Drury 1983; Blandford & Eichler 1987; 
Jones & Ellison 1991; Berezhko et al. 1996). Indeed, most shell-type 
SNRs are non-thermal radio emitters, which confirms that electrons are 
accelerated up to at least GeV energies. Moreover, the limb-brightened 
non-thermal radio emission traces the site of effective particle accelera- 
tion. 

The source SN 1006 was also the first SNR in which a non- thermal 
component of hard X-rays was detected in t he rims of the remnan t 
by ASCA iKovama et al. 1995H and ROSAT IWillingale et al. 1996H . 
whereas the interior of the remnant exhibits a thermal spectrum with 
line emission. The hard featureless power-law spectrum strongly im- 
plies a synchrotron origin of the radiation, which in turn suggests that 
electrons can be accelerated up to energies of ~ 100 TeV. Subsequent 
arcsecond resolution images by Chandra revealed a small-scale struc- 
ture in the nonthermal X-ray filame nts of the NE rim of SN 1006 
iBamba et al. 20031 , lLong et al. 20031. supp orting the idea of hig h B- 
fields in the bright limbs of the remnant I Berezhko et al. 20021 . An 
analysis of the X-ray observations from XMM-Newton by Rothenflug 
et. al (2004) leads to the conclusion that the magnetic field in the rem- 
nant is oriented in the NE-SW direction. The synchrotron emission 
would then be co ncentrated in regions where the shock is quasi-parallel 
IVolk et al. 200311 . 

Also, y-ray observations of SN 1006 were carried out by ground- 
based y-ray telescopes. A TeV y-ray signal at the level of the 
Crab flux was claim ed by the CANGAR OO-I iTanimori et al. 1998tl 
and CANGAROO-II ITan imori et al. 200111 telescopes, but subsequent 
stereoscopic observations of the source with the H.E.S.S. telescopes in 
2003 and 2004 found no evidence of VHE y-ray emission and derived 
an upper limit of $(> 0.26 TeV) < 2 4x 10~ 12 ph cm" 2 s _I at 99.9% con- 
fidence level lAharonian et al. 200511 . The CANGAROO-III telescope 
array found only an upper limit which is consistent with the H.E.S.S. 
result ITanimori et al. 2005T) . 

The initial non-detection of SN 1006 in VHE y-rays does not in- 
validate the hypothesis of nuclear particle acceleration in the shock. 
Indeed, the hadronic y-ray flux is very sensitive to the ambient gas 
density n H and hence the H.E.S.S. upper limit implies a constraint on 
na < 0.1 cnr 3 jKsenofontov et al. 2005h . Indeed, being 500 pc above 
the Galactic plane, the remnant is relatively isolated, and the gas den- 
sity around SN 1006 was recently estimated to be around 0.085 cirT 3 
[Katsuda et al. 2009]. Ksenofontov et al. (2005) furthermore showed 
that the lower limit for the VHE y-ray flux, which is given by the inverse 
Compton (IC) component derived from the integrated synchrotron flux 
and field amplification alone, was only a factor 5 below the H.E.S.S. 
upper limit. These predictions promoted deep observations with the 
H.E.S.S. telescopes. 



2. H.E.S.S. observations and analysis methods 

H.E.S.S. is an array of four 13 m diameter imaging atmospheric 
Cherenkov telescopes situated in the Khomas Highland in Namibia 
at an altitude of 1800 m above sea level IBernloehr et al. 20031 . 
iFunk et al. 2004H . The source SN 1006 was observed in 2003 with the 
two telescopes that were operational at that time and with the complete 
H.E.S.S. array in the years since. After run selection the data set com- 
prises 130 hours (live time) of observations, of which 18 hours were 
taken with two telescopes only. The latter yielded a smaller effective 
area than the data set recorded with the full array. For that reason they 
are used only in morphological studies and excluded in the spectral 
analysis. 

The data were an alysed with the Model Analysis 
Ide Naurois & Rolland 2009rl . in which shower images of all trig- 
gered telescopes are compared to a pre-calculated model by means of 
a log-likelihood minimisation. The Model Analysis does not rely on 



-42.4 



15h04m 



15h02m 



RA (hours) 



Fig. 1. H.E.S.S. y-ray significance map of SN 1006 using an in- 
tegration radius of 0.05°. The linear colour scale is in units of 
standard deviations. The white solid contours correspond to the 
regions which contain 80% of the non-thermal X-ray emission 
from the XMM-Newton flux map in the 2-4.5 keV energy range 
after smearing with the H.E.S.S. PSF, shown in the inset. The 
white dashed circles correspond to the regions that are excluded 
from background determination. 



any image-cleaning procedure and uses all pixels in the camera. The 
noise distributions in the pixels due to the night sky background are 
taken into account in the model fit and result in a superior treatment 
of shower tails. Therefore the Model Analysis results in a more 
precise reconstruction and better background suppression than more 
conventional techniques, thus leading to improved sensitivity. 

Two different sets of cuts were used: The standard cuts, including 
a minimum image charge of 60 photoelectrons (E th = 260 GeV), cover 
the full energy range and are used for the spectral analysis only. The 
hard cuts, with a larger charge cut of 200 photoelectrons, result in an 
improved signal-to-background ratio at the expense of lower statistics 
and a higher threshold of 500 GeV. These are used for the studies of 
source morphology. 

The results presented below h ave been cross-checked us- 
ing the 3D Model Analysis iLemoine-Goumard et al. 20061. 
iNaumann-Godo et al. 200911 . Both analyses yielded consistent 
results. 

Significant y-ray emission is detected from the direction of 
SN 1006, concentrated in two extended regions as shown in Fig.Q] This 
map shows the significance over a field-of-view of 1° x 1° with a pixel 
size of 0.005° obta ined with hard cu ts using the ring background sub- 
traction technique iBerge et al. 200711 and a small integration radius of 
0.05°, close to the H.E.S.S. PSF of R 6g = 0.064°. As the pixel size 
is a factor 10 smaller than the integration radius, the bins are highly 
correlated. In two regions of the map the significance of the H.E.S.S. 
observation clearly exceeds 5 <x. 

The significance distribution over the field-of-view of 2° x 2° is 
shown in Fig. [2] while Fig. [3] illustrates the area corresponding to the 
significance above a given level. The black histogram in both figures 
corresponds to the full field-of-view and exhibits strong deviation from 
a normal distribution at large significance values. The red histogram, 
restricted to the part of the field-of-view outside of the white dashed 
circles (Fig. [TJ is compatible with a normal distribution, as denoted by 
the red dashed line (Fig. O. This demonstrates that the distribution of 
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Significance 



Fig. 2. H.E.S.S. y-ray significance distribution over the full field- 
of-view ofSN 1006 (black histogram) and excluding the circular 
regions around the NE and SW emission regions ( red histogram ). 
A normal distribution (red dashed line) shows that the signifi- 
cance distribution over the rest of the field- of -view is compatible 
with expectation from statistical noise fluctuations. 
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Significance Threshold 

Fig. 3. H.E.S.S. sky area with y-ray significance above some 
threshold as a function of its value over the full fleld-of-view 
ofSN 1006 (black histogram) and excluding the circular regions 
around the NE and SW emission regions ( red histogram). 



events over the field-of-view (outside the two exclusion regions) is com- 
patible with expectation from statistical fluctuations and that systematic 
effects concerning background estimation are under control. 

3. Morphology 

Two different int egration regions were d efined a priori from the XMM- 
Newton data set IRothenfiu g et al. 200411 : a map of the flux in the 2 - 4.5 
keV energy range (to exclude thermal contamination) was smoothed 
with the H.E.S.S. PSF, and regions which contained 80% of the flux 
were calculated. The two resulting regions, denoted as NE Region and 
SW Region, are displayed as white contours in Fig.[T]and coincide well 
with the regions of largest H.E.S.S. significance. 

Excess event counts and significances for both regions are given 
in Table [T] for the two sets of cuts. The ON photons are from the re- 
gions enclosed by the solid lines in Fig. Q] while the OFF events are 
taken from regions of identical shape rotated in the field-of-view of the 
instrument around the observation position and not intersecting the ex- 
clusion regions (enclosed by dashed lines in Fig. [T}. Due to varying 



observation positions, the number of OFF regions varies from observa- 
tion to observation. Individual observation values are combined into an 
average normalisation factor (a) quoted in Table[T] Similar excess event 
counts and significances are observed in both regions, thus attesting to 
the bipolar morphology of the remnant in the TeV energy range. This 
is a highly constraining result, because due to the relatively uniform 
target density around the remnant the H.E.S.S. morphology directly re- 
flects the distribution of high-energy particles responsible for the y-ray 
emission. 



Region 


ON 


OFF 


a 


# y 


Significance 


NE, Std Cuts 


4306 


25421 


6.67 


495 


7.3 


NE, Hard Cuts 


619 


2575 


6.44 


219 


9.3 


SW, Std Cuts 


3798 


26523 


7.615 


315 


4.9 


SW, Hard Cuts 


548 


2591 


7.25 


191 


8.7 



Table 1. H.E.S.S. excess events and significances for the two re- 
gions defined from X-ray observations, a is the normalisation 
factor between OFF and ON exposures. 



Figure [4] shows the y-ray H.E.S.S. excess map, produced with 
hard cuts and the same integration radius of 0.05°, overlaid with the 
smoothed XMM-Newton flux contours. A striking similarity between 
the y-ray and X-ray emission regions is found. For a quantitative anal- 
ysis uncorrelated radial and azimuthal profiles of the H.E.S.S. excess 
events were derived and compared to the XMM-Newton profiles (Figs. 
|5]and[6]l. Again the XMM-Newton data were smoothed to match the 
H.E.S.S. point spread function, and the relative normalisation was ad- 
justed to the maximum value. Within error bars, the H.E.S.S. and 
XMM-Newton emission profiles are almost identical, thus possibly in- 
dicating a common origin. 

The geometrical X-ray centre of the SNR was derived from the un- 
smoothed XMM-Newton data by fitting them with a Gaussian radial 
profile convolved with an azimuthal profile with two Gaussian com- 
ponents, yielding 15h2m51.1s, -41d55'32.2" as the centre of the SNR 
with a radius of R = 0.239° and a thickness of dR = 0.013°. Figure [5] 
shows the radial profiles of H.E.S.S. and smoothed XMM-Newton ex- 
cess events from the centre of the SNR. When a Gaussian is fit to the 
H.E.S.S. profile (Fig. |5) the shell radius is found to be 0.24° ± 0.01° 
and the width of the radial distribution is 0.05° ± 0.01°, which is consis- 
tent with the H.E.S.S. point spread function, thereby showing that the 
emission region is compatible with a thin rim. 

The azimuthal profile, restricted to radii 0.12° < r < 0.36° from the 
centre of the SNR, is shown in Fig.|6]for H.E.S.S. data and smoothed 
XMM-Newton data in the 2 - 4.5 keV energy band. The azimuth is 
defined clockwise with zero toward the East. The H.E.S.S. profile is 
compatible with a superposition of two Gaussian emission regions al- 
most at 180° from each other, respectively centred on -143.6° ±6.1° 
(SW region) and 29.3° ± 4.0° (NE region) and with similar widths of 
33.8° ±7.0° and 27.9° ±4.0°. 

4. Spectral analysis 

Differential energy spectra of the VHE y-ray emission were derived 
for both regions above the energy threshold of ~ 260 GeV. These re- 
gions correspond to 80% of the X-ray emission (after smearing with the 
H.E.S.S. PSF) and therefore slightly underestimate the TeV emission of 
the full remnant. 

The spectra for the NE and SW regions are compatible with power 
law distributions, F(E) oc E~ r , with comparable photon indices T and 
fluxes. Confidence bands for power-law fits are shown in Fig. [7] and 
Table[2] The integral fluxes above 1 TeV correspond to less than 1% of 
the Crab flux, making SN 1006 one of the faintest known VHE sources 
(Table [2j. The derive d fluxes are well below the previously published 
H.E.S.S. upper limits I Aharonian et al. 2005ll . The observed photon in- 
dex T x 2.3 is somewhat steeper than generally expected from diffusive 
shock acceleration theory and may indicate that the upper cut-off of the 
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Fig. 6. Azimuthal profile obtained from H.E.S.S. data and XMM- 
Newton data in the 2 - 4.5 keV energy band and smoothed to 
H.E.S.S. PSF, restricted to radii 0.12° < r < 0.36° from the cen- 
tre of the SNR. Azimuth 0° corresponds to East, 90° corresponds 
to North, 180° to West and -90° to South. 



Fig. 4. H.E.S.S. y-ray image ofSN 1006. The linear colour scale 
is in units of excess counts per nx (0.05 °) 2 . Points within (0.05 ) 2 
are correlated. The white cross indicates the geometrical centre 
of the SNR obtained from XMM data as explained in the text and 
the dashed circles correspond to R± dR as derived from the fit. 
The white star shows the centre of the circle encompassing the 
whole X-ray emission as derived by Rothenflug et al. (2004) and 
the white triangle the centre derived by Cassam-Chena'i et al. 
(2008) from Ha data. The white contours correspond to a con- 
stant X-ray intensity as derived from the XMM-Newton flux map 
and smoothed to the H.E.S.S. point spread function, enclosing 
respectively 80% , 60% , 40% and 20% of the X-ray emission. 
The inset shows the H.E.S.S. PSF using an integration radius of 
0.05°. 
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Fig. 5. Radial profile around the centre of the SNR obtained from 
H.E.S.S. data and XMM-Newton data in the 2 - 4.5 keV energy 
band smoothed to H.E.S.S. PSF. 
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Fig. 7. Differential energy spectra of SN 1006 extracted from the 
two regions NE and SW as defined in Sect. [2] The shaded bands 
correspond to the range of the power-law fit, taking into account 
statistical errors. 



Region 


photon index T 


C>(> ITeV) 
(10- 12 cnrV) 


NE 
SW 


2.35±0.14. IM ±0.2 IV „ 
2.29 ± 0.1 8„ fl , ±0.2,„,„ 


0.233 ± 0.043.„«„ ± 0.047 ,„,„ 
0.155 ±0.037.„ fl , ±0.031, m 



high-energy particle distribution is being observed; however, the uncer- 
tainties on the spectrum preclude definitive conclusions on this point. 



Table 2. Fit results for power-law fits to the energy spectra. 



5. Discussion 

The source SN 1006 is an ideal example of a shell-type supernova rem- 
nant because it represents a type la supernova exploding into an ap- 
proximately uniform medium and magnetic field, thereby essentially 
maintaining the spherical geometry of a point explosion. This can be 
attributed to the fact that SN 1006 is about 500 pc above the Galactic 
plane in a relatively clean environment, where the external gas density 
is rather low, tin ~ 0.085 cnT 3 as indicated by Katsuda et al. (2009). 
Moreover, SN 1006 is one of the best-observed SNRs with a rich data- 
set of astronomical multi-wavelength information in radio, optical and 
X-rays, and all the important paramet ers like the ejected mass, it s dis- 
tance and age are fairly well-known ICassam- Chenai etal. 2008H . For 
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this reason, the semi-analytical models of Truelove & McKee (1999) 
can be approximately applied and the velocity of the shock calculated. 
The value of the shock velocity calculated by this means agrees well 
with the recent measurement in X-rays by Katsuda et al. (2009), yield- 
ing (0.48±0.04) arcsec yr~' in the synchrotron emitting regions (NE and 
SW), which corresponds to 5000 ± 400 km/s for a distance of 2.2 kpc. 
This does not contradict the value of (0.28±0.008) arcsec yr~' measured 
by Winkler et al. (2003) in the optical filaments, which are situated in 
the NW region of the remnant. All those calculations neglect the dy- 
namic role of accelerated particles however, which is potentially quite 
important. 

The basic model of VHE y-ray production requires particles accel- 
erated to multi-TeV energies and a target comprising photons and/or 
matter of sufficient density. The close correlation between X-ray and 
VHE-emission points toward particle acceleration in the strong shocks 
revealed by the Chandra observation of the X-ray filaments. Moreover, 
the bipolar morphology of the VHE-emission in the NE and SW re- 
gions of the remnant supports a major result of diffusive shock ac- 
celeration theory, according to which efficient injection of suprather- 
mal downstream charged nuclear ions is only possible for sufficiently 
small angles between the ambient magnetic field and shock normal, and 
therefore a higher density of accelerated nuclei at the poles is predicted 
HEllison et al. 199l.lMalkov & Volk 1995UV6lk et al.2003ll . 

Radio iRevnolds 1996ll and X-ray iBamba et al. 2008h data inte- 
grated over the full remnant were combined with VHE y-ray measure- 
ments to model the spectral energy distribution of the source in a simple 
one-zone stationary model. For the sake of consistency, the VHE y-ray 
energy distribution was determined from the sum of the two previously 
defined regions. In this phenomenological model the current distribu- 
tion of particles (electrons and/or protons) is prescribed with a given 
spectral shape corresponding to a power law with an exponential cutoff, 
from which emission due to synchrotron radiation, bremsstrahlung and 
IC scattering on the Cosmic Microwave Background (CMB) photons is 
computed. The n° production through interactions of protons with the 
ambient matter are obtained following Kelner et al. (2006). 

It is clear that this model oversimplifies the acceleration process 
in an expanding remnant, as discussed by e.g. Drury et al. (1989) and 
Berezhko et al. (1996). In addition one must include the uncertainties 
introduced by the dynamics of the ejecta, the nonuniform structure of 
the ambient medium and the complexities of the reaction of the accel- 
erated particles on both the magnetic field and the remnant dynamics. 
This is of importance when comparing the data to the model results 
below. 

Assuming first a purely leptonic form (Fig. [8] top), the radio and 
X-ray data constrain the synchrotron part of the SED in a way that the 
slope of the electron spectrum, which is particularly sensitive to the 
slope of the radio data, is bounded between 2.0 and 2.2, while the cutoff 
energy of electrons is limited to about 10 TeV by the X-ray data assum- 
ing a magnetic field of 30 uG. With the particle spectrum constrained 
by radio and X-ray data, the resulting magnetic field needs to be higher 
than 30 u.G so that the IC emission does not exceed the measured VHE- 
flux. A magnetic field of 30 [xG implies that assuming Bohm diffusion, 
electrons of 1 TeV are confined in a shell of the width of 10 arcseconds, 
which is much smaller than the PSF of the H.E.S.S. instrument and is 
therefore compatible with the radial profile shown in Fig. [5] However, 
while this simple leptonic scenario can account for the measured VHE 
y-ray flux, it fails to reproduce the slope of the VHE spectrum, which is 
much harder than the expectations from the IC process (see Fig.[8]top). 
But it should be noted that non-linear Fermi shock acceleration as re- 
viewed by Malkov & Drury (2001) usually predicts curved cosmic ray 
spectra with different spectral shapes for protons and electrons. There is 
a hint of spectral curvature observed in t he case of Tycho's and Kep ler's 
supernova remnants in the radio regime (Reynolds & E llison 19921] ■ F° r 
SN 1006 there is also an indication of th e curvature of the e lectron spec- 
trum in the GeV to TeV energy range lAllenetal. 2008H . These non- 
linear effects, which also might well introduce a spectral curvature in 
the VHE regime, are not addressed by this simple model. 

In a second dominantly hadronic model (Fig. [8] middle) TeV emis- 
sion results from proton-proton interactions with 7r° -production and 
subsequent decay, whereas the X-ray emission is still produced by lep- 
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Fig. 8. Broadband SED models of SN 1006 for a leptonic 
scenario (top), a hadronic one (centre) and a mixed lep- 
tonic/hadronic scenario (bottom). Top: Modelling was done by 
using an electron spectrum in the form of a power law with an 
index of 2.1, an exponential cutoff at 10 TeV and a total en- 
ergy ofW e — 3.3 X 10 47 erg. The magnetic field amounts to 
30 /uG. Centre: Modelling using a proton spectrum in the form 
of a power law with an index of 2.0, an exponential cutoff at 
80 TeV and a total proton energy o/W p = 3.0 X 10 50 erg (us- 
ing a lower energy cut off of 1 GeV). The electron/proton ratio 
above 1 GeV was K ep = 1 X 10~ 4 with an electron spectral index 
of 2.1 and cutoff energy at 5 TeV. The magnetic field amounts to 
120 fiG and the average medium density is 0.085 cmT^ . Bottom: 
Modelling using a mixture of the above two cases. The total pro- 
ton energy was W p = 2.0 x 10 50 erg, with K ep = 7x 10~ 3 , with 
exponential cutoffs at 8 TeV and 100 TeV for electrons and pro- 
tons resp ectively. The ma gnetic field a mounts to 45 [iG. Th e ra- 
dio data /[Reynolds 199&1 . X-ray data kBamba et al. 2008\1 and 
H.E.S.S. data (sum of the two regions) are indicated. The follow- 
ing processes have been taken into account: synchrotron radia- 
tion from primary electrons (dashed black lines), IC scattering 
(dotted red lines), bremsstrahlung (dot-dashed green lines) and 
proton-proton interactions (dotted blue lines). The Fermi/LAT 
sensitivity for one year is shown (pink) for Galactic ( upper) and 
extragalactic (lower) background. The latter is more representa- 
tive given that SN 1006 is 14° north of the Galactic plane. 
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tonic interactions. A rough representation of the effect of spectral cur- 
vature is included by allowing for a slightly harder spectral index for 
protons than for radio-emitting electrons. A lower electron fraction al- 
lows us to account for the X-ray and radio emission with a higher field 
value of 120 (xG, which is consistent with magnetic field amplification 
at the shock, as indicated by the above-mentioned measurements of thin 
X-ray filaments. Assuming an average medium density of 0.085 cm -3 
and a proton spectral index of 2.0 with a cutoff energy of 80 TeV (in- 
ferred from the maximum energy of TeV photons), this model requires 
a high overall fraction of about 20%, of the supernova energy to be 
converted into high-energy protons. Here E$ N = 1.4 x 10 51 erg was as- 
sumed, near the upper end of the typical range of type la SN explosion 
energies (e.g. Woosley et al. 2007), as the assumed density, observed 
radius and known age of SN 1006 appear to require a higher than aver- 
age explosion energy. Given that the VHE emission is concentrated in 
polar regions of the shell, the local shock acceleration efficiency would 
then be several times higher than this fraction. 

In a third example (mixed model), hadronic and leptonic processes 
contribute almost equally to the very high-energy emission. The elec- 
tron spectrum is similar to the aforementioned leptonic case and the 
total proton energy is set to 14% of the mechanical supernova energy 
with the electron/proton ratio K cp = 3.9 x 10~ 3 , thus leaving the mag- 
netic field and the cutoff energy of protons the only free parameters. In 
the example shown in Fig.[8](bottom panel) the magnetic field amounts 
to 45 uG and the cutoff energy of protons is 100 TeV. This example 
illustrates that in this simple one-zone case it is possible to reproduce 
all the multi- wavelength data on SN 1006 to a reasonable degree of 
accuracy including the slope of the VHE-data. While these considera- 
tions cannot exclude any of the astrophysical scenarios, they serve as a 
quantitative illustration of the various alternatives. 

Values of total electron and proton energy, cutoff energy and mag- 
netic field obtained in the three aforementioned cases are summarised 
in Table[3] These parameters yield very similar values when the NE and 
SW regions are adjusted independently. 
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10 




3.3 




30 


Hadronic 
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80 


0.3 


3.0 


120 


Mixed 


8 


100 


1.4 


2.0 


45 



Table 3. Parameters used in the spectral energy modelling 
shown in Fig. [3] Spectral indices have been fixed to 2.1 and 2.0 
respectively for electrons and protons. 



M ore elaborate models using e.g. a nonlinear kinetic acceleration 
theory jBerezhko et al. 20"09tl go beyond the simple approach developed 
here and lead to precise predictions that could be quantitatively tested 
against the data. Several effects which were not included in the simple 
model above would alter the total energy in accelerated particles re- 
quired for the hadronic component. Beyond the spectral curvature men- 
tioned previously, these include the higher compression of the target 
matter induced by the dynamical reaction of the accelerated particles, 
and consideration of the heavier nuclei composition of the accelerated 
hadrons instead of the pure protons assumed here. Measurements in the 
GeV-energy range would be pivotal to distinguish between the differ- 
ent scenarios. Unfortunately, the sensitivity of the Fermi Large Area 
Telescope for one year as given in Atwood et al. (2009) is of a factor 
of the order of 10 too low (depending on the model and the exact dif- 
fuse background flux) to measure the predicted flux at 1 GeV as shown 
in Fig(8] which makes the detection of SN 1006 by Fermi LAT rather 
unlikely. 

6. Conclusions 

Very high energy y-rays from SN 1006 have been detected by H.E.S.S. 
The measured flux above 1 TeV is of the order of 1% of that detected 



from the Crab Neb ula and therefore comp atible with the previously pub- 
lished upper limit | Aharonian et al. 200511 . The bipolar morphology ap- 
parent in y-rays is consistent with the non-thermal emission regions also 
visible in X-rays. As the VHE-shell is compatible with a scenario of thin 
rim emission, particle acceleration in the very narrow X-ray filaments, 
which are signatures of shocks, is also likely to be at the origin of the 
y-ray signal. The measured flux level can be accounted for by inverse 
Compton emission assuming a magnetic field of about 30 /jG. A mixed 
scenario including leptonic and hadronic processes and taking into ac- 
count the ambient matter density estimated from observation also leads 
to a satisfactory description of the multi-wavelength spectrum, assum- 
ing a high proton-acceleration efficiency. None of the models can be 
excluded at the level of modelling presented here. 
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